Stochastic resonance in a model of opinion formation on small-world networks 
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We analyze the phenomenon of stochastic resonance in an Ising-hke system on a small- 
world network. The system, which is subject to the combined action of noise and an external 
modulation, can be interpreted as a stylized model of opinion formation by imitation under 
the effects of a "fashion wave" . Both the amplitude threshold for the detection of the external 
modulation and the width of the stochastic-resonance peak show considerable variation as 
the randomness of the underlying small-world network is changed. 

PACS numbers: 87.23.Ge, 89.65.-s, 05.40.-a 
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INTRODUCTION 

Many mathematical models of social processes, in- 
spired by analogies with physical systems, have recently 
been formulated in order to describe a wide spectrum of 
phenomena 0|. Economic and financial processes, dis- 
ease spreading and information propagation, and evolu- 
tion of social structures, among others, have been ana- 
lyzed along such lines. Weidlich has proposed a model 
of public opinion formation which, in the simplest vari- 
ant, considers that the individuals can adopt two differ- 
ent opinions. The opinion of a given individual may be 
influenced by that of the neighbors, making it to change 
with a certain probability. In the present work, we ana- 
lyze a similar model of social influence and opinion for- 
mation. In contrast with previous studies of this class 
of models, however, we are interested at taking into ac- 
count the network of social interactions. Specifically, we 
model the underlying social structure of the population 
as a small-world network. Small-world networks in- 
corporate two main features of real social interactions. 
First, they are highly clustered, which means that any 
two neighbors of a given site have a relatively large prob- 
ability of being in turn mutual neighbors. Second, the 
mean number of intermediaries between any two sites is 
quite small, and increases very slowly as the total num- 
ber of sites grows [||. This is precisely the small- world 
property, originally discussed by Milgram as a typical fea- 
ture of social communities and relationships Small- 
world networks can be considered as partially disordered 
structures. The construction procedure makes it possi- 
ble to control the degree of disorder of a given network, 
ranging from ordered lattices to completely disordered 
graphs. Among many other applications, these networks 
have been used as a model of social structures in the 
study of disease |§, ||, |[ 0, ^ and rumor propagation [|| . 

As pointed out in the next section, our model is an 
Ising-like system whose dynamics is driven by the ma- 
jority rule. Equilibrium properties of Ising systems on 
small-world networks have been studied in detail. It has 
been shown that, in asymptotically large systems, even 



an infinitesimal amount of disorder on a one-dimensional 
Ising lattice is sufficient for the system to undergo a ferro- 
magnetic phase transition as the temperature varies 
Here, instead, we consider the case where the system is 
maintained out of equilibrium by an external modula- 
tion. In the frame of the problem of opinion formation, 
this modulation plays the role of a periodic "fashion" 
wave, as actually observed in certain real situations 
The combination of this modulation with noise — which, 
in equilibrium, would play the role of the temperature — 
leads naturally to the consideration of stochastic reso- 
nance |11 . In the phenomenon of stochastic resonance, 
an enhancement of the response of a system to an ex- 
ternal modulation is obtained by an adequate choice of 
the noise level. This effect has attracted great interest 
due to both its potential technological applications and 
its connection with biological detection mechanisms. The 
range of possible applications in science includes paleocli- 
matology Jl^ , electronic circuits p^, lasers jl^, ^ , 
chemical |17 and biological systems [|18|, [l9| . Stochas- 
tic resonance in Weidlich's model of opinion formation 
has already been studied In the present work, we 

analyze this phenomenon in a closely related model of 
social imitation, incorporating a small-world network as 
the underlying structure of social interactions. 



THE MODEL 

As advanced in the Introduction, we consider an Ising- 
like model on a network with a controlled degree of dis- 
order, namely, a small-world network. The network is 
constructed as follows [^. We start with an A^-node 
one-dimensional ordered lattice with periodic boundary 
conditions, where each node is linked to its 2K nearest 
neighbors {K clockwise and K counterclockwise). The 
network is then disordered, rewiring each of the K clock- 
wise connections of each node z to a randomly chosen 
node j , with probability p. Double and multiple connec- 
tions are forbidden, and realizations where the network 
becomes disconnected are discarded. The final degree of 
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FIG. 2: Response / as a function of the modulation amplitude 
A, for rj = 0.1, and four values of the randomness p. 



FIG. 1: Normalized histograms of the residence time tr, for 
p = 0.8. Three different noise intensities rj are considered. 
Each histogram has been constructed recording the residence 
times over an interval of 8 x 10^ time units, with a subsequent 
average over 20 realizations of the small-world network and 
the initial condition. The insert in the third plot shows the 
same data with logarithmic scale in the vertical axis. 

disorder, or randomness, is determined by the probability 
p. As a result of the disordering process, some shortcuts 
between otherwise distant regions are created. Note that, 
independently of the value of p, the average number of 
hnks per site is always 2K. 

Each site of the small-world network is occupied by an 
individual with two possible states, fii = ±1. At each 
time step, an individual i is selected at random from 
the whole population, and the following three evolution 
rules are successively applied to its state fii. First, fii is 
changed according to the majority rule, i.e. 

fii ^ sign^ fij , (1) 

where the sum runs over the neighborhood of site i. If 
the sum equals zero, /i^ is left unchanged. Second, we 
apply the external modulation. The effect of this mod- 
ulation is prescribed by a harmonic function of time, 
/(t) = ^cos(27rt/T), with < A < 1. Namely, the state 
of individual i changes according to 

fii sign f{t), (2) 

with probability Ilr{t) = |/(^)|- Finally, noise is applied 
as follows. With constant probability n„ the state 
is left unchanged. With the complementary probability 
1 — n„, fii is assigned one of the two possible values ±1, 
chosen at random with probability 1/2. We call = 
1 — n„ the noise intensity. 



In the absence of modulation, at low noise levels, and 
for moderate and large randomness p, our population 
evolves towards an almost homogeneous state, where fii 
is the same at practically all sites. The average of fii 
over the population, Jl = '^^^-^ fii, is close to ±1. 
Due to the effect of noise, occasional transitions between 
states with opposite signs are observed. For sufficiently 
small p, on the other hand, Jl fluctuates around zero 
with relatively small amplitude. These two regimes cor- 
respond, respectively, to the ferromagnetic and paramag- 
netic phases predicted for equilibrium small-world Ising 
systems In our simulations, we will mainly consider 
randomness and noise levels for which the nonmodulated 
system is in the ferromagnetic-like phase. All the cases 
with p > will in fact correspond to that phase. In 
this situation, the noiseless version of our system exhibits 
bistability. The paramagnetic-like phase will only be en- 
countered in the special case p = 0. 

The effect of the external modulation on the 
ferromagnetic-like phase consists of a periodic modifica- 
tion in the transition probability between the two states 
with /I ~ ±1. In the paramagnetic- like phase, on the 
other hand, Jl oscillates around zero at the same fre- 
quency as the modulation. In both cases, of course, the 
modulation is combined with the random fluctuations of 
noise. It is this nontrivial interplay, which may result in 
the enhancement of the system response, that we focus 
on in our numerical simulations. 

NUMERICAL RESULTS 

We have performed extensive numerical simulations of 
the model presented in the previous section, considering 
small-world networks of iV = 10^ nodes with K = 3. The 
time interval assigned to an evolution step is At = 
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FIG. 3: The response /, in logarithmic scale, as a function of 
the inverse modulation amplitude, , for four values of p. 
Straight lines are least-square fittings to the linear regime of 
each data set. 

SO that, on the average, the state of each site is updated 
once during each time unit. Each realization starts with 
the generation of the random network and the random 
initialization of the state of all the elements. The evolu- 
tion is analyzed after a transient of 10^ time units, during 
which the system is found to reach a steady large-time 
behavior. In particular, we record the residence times 
tr, i. e. the intervals between two consecutive changes 
of sign in the average state JI. Note that the external 
modulation is expected to induce changes in the sign of 
JI with a typical period close to r/2. The values of tr are 
used to construct a histogram, showing the frequency of 
occurrence of each residence time. Throughout all the 
calculations we considered tr — 400. 

Figure ^ illustrates this histogram for fixed modulation 
and small-world randomness, and three levels of noise. 
For low noise, the histogram shows a peak at tr — t/2, 
as expected, but several odd harmonics, tr = A;r/2 with 

= 3, 5, 7, . . ., are also observed. These harmonics corre- 
spond to residence intervals during which the modulation 
fails, one or more times, at inducing the transition. At 
intermediate noise levels, corresponding to the resonance 
intensity, only the pick around tr — t/2 is important. 
For higher noise levels the transitions are disordered and 
highly frequent, so that the residence times are usually 
very short and the histogram shows a peak near tr — 0. 
We stress that the analysis of these patterns has been 
employed as a quantitative description of stochastic res- 
onance in other systems ^ , as an alternative to the 
study of the signal-to- noise ratio In fact, a measure 
of the sensitivity to the external modulation is given by 
the area / under the peak around t/2 in the (normalized) 
histogram. Below, we refer to the area / as the response 
of the system. 
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FIG. 4: The response / as a function of the noise intensity 
r], for three values of the small-world randomness p. Lines 
joining dots are drawn as a guide to the eye. 

First, we analyze the behavior of our system as the 
modulation amplitude A is varied. Fixing the intensity 
of noise rj we observe, for p > 0, the existence of a thresh- 
old in A below which the system is not able to respond 
to the modulation. We find that, as the randomness p 
increases, the threshold amplitude grows as well. Figure 
H shows the response / for a noise intensity = 0.1 as 
a function of A, for several values of p. The nature of 
this threshold phenomenon is clarified in Fig. ^, where 
we plot / in logarithmic scale as a function of the inverse 
amplitude for three values of p. The approximately linear 
dependence for small / suggests the functional form 

/(A) =/oexp(-AoM). (3) 

This form is reminiscent of Kramers formula for the tran- 
sition frequency lo between equilibrium states under the 
action of thermal noise T, lo (x exp(— Tq/T), where To 
is related to the potential barrier separating the equilib- 
ria In our case, for a fixed noise level, the role of 
temperature — as the factor that induces the transition — 
is played by the external modulation. Exploiting this 
analogy, our model can be thought of, at the macro- 
scopic level, as a system subject to an effective bistable 
potential, with a barrier proportional to the slope Aq . It 
is apparent from Fig. || that the slope depends on the 
small-world randomness. As p grows from 0.1 to 0.9, the 
slope increases by a factor of about 2.5, indicating that 
the minima in the effective potential become deeper for 
larger randomness. 

Taking into account the results shown in Figs. |^ and 
^, in our study of dependence of the response / on the 
noise intensity rj we fix A — 0.2 which, for the values of p 
and 77 considered in our simulations, is always above the 
detection threshold. In Fig. ^, we show / as a function of 
eta for several values of p. The existence of a resonance 
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FIG. 5: Area S under the stochastic-resonance curves (Fig. 
^ and lengths Lm for several values of M (see main text), as 
a function of the randomness p. Lines joining dots are drawn 
as a guide to the eye. 



as the noise intensity is varied is apparent. At the same 
time, an effect related to the variation of the randomness 
p is observed. Namely, the higher is p, the broadest is 
the resonant interval in the eta. In other words, accurate 
tuning of the noise level in order to enhance the response 
is less important in more disordered networks. 

Figure ^ presents a quantitative characterization of the 
broadening of the resonance peaks as p grows. First, we 
show the area S of the peaks as a function of p. Moreover, 
we show the length Lm of the intervals of noise intensity 
where the response / is above the level Im = O.OIM, for 
several values of M. 



CONCLUSION 

We have presented an Ising-like model of opinion for- 
mation by imitation on a small-world network, subject 
to the action of an external modulation and noise. The 
study is focused on the phenomenon of stochastic reso- 
nance, which underlies the macroscopic behavior of this 
system. The main control parameters in our analysis are 
the amplitude A of the external modulation, the intensity 
of noise rj, and the small-world randomness p. The first 
aspect studied here is the existence of a threshold in the 
amplitude of the external signal, below which the system 
is not able to detect the modulation. This threshold is 
higher as the disorder of the network increases, and dis- 
appears for p ~ 0. This fact suggests that the existence 
of a threshold in A is directly related to the bistable na- 
ture of the effective potential which acts on the system 
in its ferromagnetic-like phase, i. e. for moderate and 
large randomness. It is interesting to mention that the 
increase of the threshold amplitude with the randomness. 



which amounts to a growth in the effective potential bar- 
rier, has a well-known correlate in social imitation pro- 
cesses subject to the action of noise |Q. The resistance 
to collective changes of opinion grows drastically when 
long-range social interactions are allowed. A population 
restricted to local interactions, instead, is susceptible to 
global opinion transitions by diffusive-like propagation. 

As expected, the analysis of the response of the system 
as a function of the noise intensity rj — provided that the 
amplitude A is above the detection threshold — reveals 
the occurrence of a resonance phenomenon. Though the 
response of the system at the resonance is practically in- 
dependent of the randomness p, an enhancement of the 
phenomenon is apparent from the broadening of the res- 
onance curves as a function of rj (Figs. ^ and ||). The 
convenience of broadening the stochastic-resonance peak 
to avoid noise tuning has been discussed in connection 
with a variety of systems . 

The existence of stochastic resonance in a model of 
opinion formation yields the appealing implication that 
there is an optimal noise level for a population to respond 
to an external "fashion" modulation . Lower noise in- 
tensities lead to the dominance of the majority's opinion, 
irrespectively of the external influences, while sufhciently 
stronger fluctuations prevent the formation of a definite 
collective opinion. We have here shown that, in this phe- 
nomenon, the underlying structure of social interactions 
plays a role, allowing for a looser tuning of the resonance 
noise intensity as the disorder grows. 

Let us finally mention that we have tested the robust- 
ness of these results by studying a few variations of our 
model. For instance, we have considered the case where 
the external modulation acts only when it is not possible 
to define an individual's state by the majority rule, i. e. 
when the sum in Eq. (Q) equals zero. In this case, the in- 
dividual adopts the state dictated by the "fashion" wave 
with probability Ilr{t), defined in Sect. II, or preserves 
the previous state with the complementary probability. 
In all cases we got consistent results, qualitatively similar 
to those presented here. Though our model is an oversim- 
plified caricature of any real process, the results illustrate 
the complex interplay between fluctuations, external in- 
fluences, and interaction structure that is expected to 
take place in actual societies. 

The authors thank H. S. Wio for fruitful suggestions. 
MK thanks Fundacion Antorchas for financial support. 



kupermanQcab.cnea.gov .ai 



zanetteQcab. cnea.gov. ai 



Electronic address: 
Electronic address: 
W. Weidlich. Phys. Rep. 204, 1 (1991). 
D. J. Watts and S. H. Strogatz, Nature 393, 440 (1998). 
A. Barrat and M. Weigt, Eur. Phys. J. B 13 (2000) 547- 
560. 

S. Milgram, Psychol. Today 2, 60 (1967). 



5 



[5] M. Kupcrman and G. Abramson, Phys. Rev. Lett. 56, 
1906 (2001). 

[6] C. Moore and M. E. J. Newman, Phys. Rev. E 60 (1999) 
5678. 

[7] R. Pastor-Satorras and A. Vespignani, Phys. Rev. E 63, 
066117 (2001). 

[8] N. Zekri and J. P. Clerc, Phys. Rev. E 64, 056115 (2001). 
[9] D. H. Zancttc, Phys. Rev. E 64 050901(R) (2001). 
[10] P. Babincc, Phys. Lett. A, 225, 179 (1997). 
[11] L. Gammaitoni, P. Hanggi, P. Jung and F. Marchesoni, 

Rev. Mod. Phys. 70 223 (1988). 
[12] R. Benzi, A. Sutera and A. Vulpiani, J. Phys. A 14, L453 
(1981). 

[13] S. Fauve and F. Heslot, Phys. Lett. A 97, 5 (1983). 
[14] R. N. Mantegna and B. Spagnolo, Phys. Rev. E 49, 

R1792 (1994). 

[15] J.M. lannelh, A. Yariv, T.R. Chen and Y.H. Zhuang, 
Appl. Phys. Lett. 65, 1983 (1994). 



[16] A. Simon and A. Libchaber, Phys. Rev. Lett. 69, 3375 
(1992). 

[17] V. Potrov, Q. Ouyang and H. L. Swinney, Nature 388, 
655 (1997). 

[18] J. K. Douglas, L. Wilkens, E. Pantazelou and F. Moss, 
Nature 365, 337 (1993). 

[19] S. M. Bezrukov and L Vodyanoy, Nature 378, 362 (1995). 

[20] F. Castro, M. Kuperman, M. Fuentes and H. Wio, Phys. 
Rev. E 64, 0511058 (2001). 

[21] G. W. Gardiner, Handbook of Stochastic Methods for 
Physics, Chemistry and the Natural sciences (Berlin, 
Springer, 1985). 

[22] H. Peyton Young, Individual Strategy and Social Struc- 
ture (Princepton University Press, Princeton, 1998). 

[23] S. Bouzat and H. S. Wio, Phys. Rev. E 59, 5142 (1999). 

[24] G. J. Tessone, H.S. Wio and P.Hanggi, Phys. Rev. E 62, 
4623 (2000). 



